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SECTION  I:  PLASMA  THEORY  &  SIMULATION 

A.  Electrostatic  Ion-Ion  Two  Streaming  Instability 
in  a  Thermal  Barrier  Cell 

V.  A.  Thomas  (Dr.  W.  M.  Nevins,  LLNL) 

A  report  on  this  project  is  in  preparation  for  publication. 


B.  Low  Frequency  Electromagnetic  Ion-Streaming  Instabilities 

V.  A.  Thomas 


I.  Introduction 

Many  studies  have  been  directed  at  instabilities  arising  from  ion  beams  counter  streaming 
along  the  field  lines  as  in  Refs.  1*2  and  the  references  in  Ref.  2.  In  particular,  low  frequency 
electromagnetic  instabilities  may  arise  in  the  following  situations;  (1)  neutral  beam  injection 
parallel  to  the  magnetic  field,  (2)  runaway  ions  in  a  tokamak  as  in  Ref.  3,  (3)  steady-state 
tokamaks  sustained  with  neutral  beam  injection  as  in  Ref.  4,  (4)  the  counter  streaming  torus 
reactor  concept  as  in  Ref.  S,  and  (S)  various  space  physics  applications  such  as  shocks.  Ion- 
beam  instabilites  may  also  be  important  in  some  tandem  mirror  designs.  Simulation  studies  of 
low  frequency  electrostatic  streaming  instabilities  are  quite  numerous  whereas  the  same  cannot 
be  said  of  low  frequency  electomagnetic  streaming  instabilities.  Our  goal  is  to  examine  some 
aspects  of  some  of  the  situations  indicated  above. 

II.  Preliminary  Results 

In  this  section  we  present  some  initial  results  for  the  restricted  case  of  propagation  parallel 
to  the  ambient  magnetic  field.  The  simulation  code  has  fully  nonlinear  ions  and  fluid  electrons 
and  uses  the  Darwin  approximation  together  with  quasineutrality  to  solve  for  the  field  quantites. 
The  code  was  developed  by  N.  Otani  from  an  algorithm  from  D.  Hamed  as  in  Ref.  6.  The  code 
is  able  to  solve  2-Va  dimensional  problems  with  the  ambient  magnetic  field  in  an  arbitrary  direc¬ 
tion.  For  these  test  cases  we  use  the  code  in  a  l-d,3  v  periodic  form. 

Our  test  case  consists  of  the  problem  considered  in  Ref.  7,  the  problem  of  cold  counter- 
streaming  ion  beams  with  arbitrary  relative  density.  The  streaming  direction  is  the  same  as  the 
direction  of  the  ambient  magnetic  field  and  we  use  the  center*of-momentum  frame  of  refer¬ 
ence.  The  plasma  is  also  charge  neutral.  Then  making  the  Darwin  approximation  of  negligible 
transverse  displacement  current  and  neglecting  the  electron  inertial  effects  we  obtain,  as  in  Ref. 
7 
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Here  c  is  the  speed  of  light,  tj  is  the  ratio  of  the  less  dense  ion  beam  particle  density  to  the 
electron  density,  is  the  ion  plasma  frequency  calculated  using  the  total  ion  density  ,  and  U 
is  the  relative  drift  between  the  two  ion  beams. 

In  Fig.  1  and  Fig.  2  we  compare  linear  growth  rates  and  real  frequencies  for  two  cases 
from  Ref.  7.  The  agreement  between  simulation  and  theory  is  quite  good.  The  parameters  for 
our  simulations  were  flaAt  -  0.01  ,  ng  -  64,  and  np  »  2000  .  The  beams  were  initially  cold 
and  only  the  first  Fourier  mode  was  retained. 


III.  Future  Directions 

In  the  future  we  plan  to  examine  the  nonlinear  consequences  of  some  of  these  instabili¬ 
ties  and  the  various  processes  important  for  stability.  Also  it  may  be  necessary  to  perform 
simulations  that  include  propagation  at  an  angle  to  the  magnetic  field  so  as  to  the  include  the 
instability  with  the  lowest  threshold  as  explained  in  Ref.  1.  This  unstable  mode  is  an  Alfven 
mode  in  contrast  to  the  mode  presented  in  this  paper  which  is  a  magnetosonic  mode.  The 
saturation  levels  for  these  two  modes  need  to  be  examined  carefully. 
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SECT10N  I:  PLASMA  THEORY  &  SIMULATION 


C.  Aifven  Ion  Cyclotron  Instability  Particle  Simulations 

Niels  F.  Otani  (Dr.  G.  R.  Smith,  LLNL) 


Details  of  the  progress  made  in  our  study  of  the  Alfvln  ion-cyclotron  instability  for  the  period 
covered  by  this  publication  will  be  presented  as  part  of  a  later  report,  currently  in  preparation. 
Four  constants  are  calculated  for  the  motion  of  ions  in  a  single,  purely  transverse  wave  propagating 
along  a  uniform  magnetic  held.  A  simple,  physically  intuitive  derivation  of  the  A1C  dispersion 
relation  is  presented  using  an  infinite  fluid  representation  of  the  standard  Vlasov  phase-space  fluid. 
These  results  are  used  in  the  analysis  of  computer  simulations  of  the  instability.  Simulations  to 
this  point  have  employed  "four-spokes"  ion  distribution  functions  similar  to  those  described  in  a 
previous  QPR,1  generalized  by  the  inclusion  of  several  parallel  velocity  beams  to  model  finite  Tq\. 
Two  separate  simulation  codes  are  being  used.  TRACY,  a  one-dimensional  magneto-inductive  code 
has  already  been  described.1  Additionally,  a  new,  two-dimensional  quasineutral  Darwin  algorithm 
has  been  written,  and  will  be  detailed  in  the  next  QPR.  Currently,  this  code  is  being  used  to 
check  the  results  being  produced  by  TRACY,  and  is  therefore  being  run  as  a  one-dimensional 
code.  Linear  growth  rates  observed  in  both  codes  agree  well  with  those  predicted  theoretically 
for  continuous  bimaxweliian  distributions  for  %  1,  in  the  vicinity  of  the  peak  of  the  growth 
spectrum.  Discrepancies  in  the  growth  rates  for  larger  values  of  k  and  smaller  0i±  are  discussed  with 
the  help  of  numerical  solutions  to  the  dispersion  relation  for  I)  multibeam  distribution  functions 
and  2)  continuous  bimaxweliian  distribution  functions.  Observations  of  particle  trapping  and  other 
modifications  to  the  ion  distribution  near  wave  saturation  are  explained  in  terms  of  the  approximate 
constants  of  the  particle  motion. 

1  Third  Quarter  Progreee  Report  on  Platma  Theory  and  Simulation.  July  1,  1081  -  Sept.  30, 
1981. 


SECTION  I:  PLASMA  THEORY  &  SIMULATION 


D.  Theory  of  Asymmetric  Double  Layers 

K.  Y.  Kim 

An  ERL  Report  Memo  with  this  title  exists,  with  abstract  as  follows: 

We  present  analytic  solutions  for  asymmetric  double  layers  which  satisfy  the  time  stationary 
Vlasov-Poisson  system  and  which  require  a  double  valued  Sagdeev  potential  as  a  function  of 
physical  potential:  it  is  pointed  out  that  any  distribution  function  having  an  analytic  density 
representation  as  a  polynomial  power  series  of  potential  can  never  satisfy  the  asymmetric  dou¬ 
ble  layer  boundary  conditions.  K-dV  like  equation  is  found  showing  a  relationship  between  the 
speed  of  asymmetric  double  layer  and  the  degree  of  asymmetry. 


£.  Kinetic  Equations  and  the  Effects  of  Electron-PIasmon  Coupling 
on  Electromagnetic  Wave-Plasma  Interaction 

K.  Y.  Kim  (Prof.  C.K.  BirdsaU  and  Dr.  T.L.  Crystal) 

An  ERL  Report  with  this  title  is  in  preparation,  with  abstract  as  follows: 

We  have  obtained  the  coupled  equations  for  the  electron  and  plasmon  distribution  func¬ 
tions,  taking  into  account  the  multiphoton  effect  and  the  nonlinear  enhancement  of 
plasmons  in  the  presence  of  an  electromagnetic  wave.  By  using  a  time  dependent  unitary 
transformation,  the  nonlinear  damping  rate  by  electron-plasmon  coupling  is  calculated 
non-perturbatively  with  respect  to  the  electromagnetic  field  intensity. 
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F.  Plasma  Surface  Interactions 

Prof.  C.  K.  Birdsall 

(work  done  at  IPP,  Nagoya,  Japan  and  Cal  Tech,  Pasadena) 

(1)  Introduction:  Plasmas  and  surfaces  (walls)  interact  in  a  great  variety  of  ways.  The 
plasma  near  the  wall  forms  a  sheath  which  may  be  rather  far  from  electrically  neutral  relative  to 
the  body  of  the  plasma.  The  sheath  region  may  extend  a  few  Debye  lengths  into  the  plasma,  or 
a  few  hundred  Debye  lengths.  The  wall  may  absorb  and  emit  particles  in  a  variety  of  ways  and 
may  also  be  coupled  to  external  sources  and  circuits,  both  active  and  passive.  The  particle 
penetration  in  the  wall  may  range  from  zero  to  a  few  atomic  layers  to  many  layers.  Knowledge 
of  the  physics  and  chemistry  of  the  volume  and  surface  interactions  is  of  interest  scientifically 
and  for  applications  such  as  fusion  first  walls,  space  ship  charging,  and  industrial  processing. 

The  current  state  of  knowledge  of  plasma-surface  interactions  may  be  termed  fair  theoret¬ 
ically  and  experimentally  with  a  very  extensive  literature,  but  yet  not  wholly  satisfactory.  What 
happens  to  the  standard  models  with  Te »  T,  when  the  ions  are  made  warm,  with  T,  compar¬ 
able  to  T„  and/or  KT  is  in  the  keV  ranges  as  in  fusion  devices?  What  is  the  detailed  volume 
and  surface  physics  and  chemistry  in  plasma  ("dry”)  etching?  Are  the  time-independent  and 
low-frequency  answers  sufficiently  close  to  the  time  average  of  the  time-dependent  behavior, 
especially  in  light  of  known  large  amplitude  nonlinear  and  possibly  turbulent  sheath  behavior? 
There  are  many  questions  yet  to  be  answered. 

Computer  simulation  using  particles  (and/or  fluids  also,  where  appropriate)  is  well  suited 
to  explore  the  time-dependent  plasma  behavior  in  plasma-wall  interactions,  where  the  plasma  is 
not  neutral,  non-Maxwellian,  possibly  in  a  large  amplitude  state,  with  both  fluid-like  and 
particle-like  waves  existing  together.  The  wall  behavior  may  also  be  explored  using  particles  for 
example;  the  plasma-surface  model  may  be  considered  as  two  or  more  adjacent  plasmas,  each 
with  characteristic  scale  lengths  and  times,  having  electrons  and  ions  and  neutrals  (atoms, 
molecules). 

(2)  Acknowledgements:  This  research  is  relatively  new  to  us  and  was  begun  in  coopera¬ 
tion  with  the  Institute  of  Plasma  Physics,  Nagoya  University,  Nagoya,  Japan.  C.  K.  Birdsall  was 
at  IPP  from  September  1,  1981  until  February  28,  1982  (on  leave),  with  a  short  visit  in  Sep¬ 
tember  1982,  working  with  Dr.  T.  Kamimura  and  programmer  Y.  Ohara  on  plasma  sheath 
simulations  during  the  second  half  of  this  period,  as  reported  herein.  Prof.  Birdsall  is  grateful 
to  Dr.  Kamimura  for  being  a  gracious  host  and  for  providing  the  excellent  services  of  Mr. 
Ohara,  who  made  setting  up  the  sheath  simulation  work  a  real  pleasure,  and  for  partial  support 
from  the  Japan  Ministry  of  Education,  through  IPP.  We  plan  on  continuing  joint  research  in 
this  area  of  plasma-wall  interactions,  as  a  long-term  cooperative  project,  with  additional  person¬ 
nel  exchanges  and  including  the  possibility  of  a  computer  link,  Nagoya  to  California.  Part  of 
the  diagnostics  of  the  results  was  done  at  California  Institute  of  Technology,  with  thanks  to  my 
hosts  Profs.  W.  B.  Bridges,  R.  W.  Gould,  and  to  the  Chevron  Research  Co.  for  my  support  as 
Visiting  Professor,  Brian  Leahy  did  the  programming. 

(3)  Related  Areas  -  Bibliography:  Plasma-surface  interactions  with  large  and  fluctuating 
potentials,  tend  to  be  related  to  other  examples  of  large  potentials  in  plasmas,  such  as:  forma¬ 
tion  and  fluctuations  in  double  layers  and  in  virtual  cathodes  and  virtual  anodes;  free  expan¬ 
sion  of  a  warm  plasma;  launching  and  propagation  and  decay  of  very  large  amplitude  waves 
(BGK,  solitons,  etc.);  thermal  barriers  (really  potential  barriers  or  wells  in  tandem  mirror 
fusion  devices).  The  author  took  advantage  of  the  excellent  library  and  library  personnel  at 
IPP,  Nagoya  to  locate  a  representative  group  of  references  in  these  related  areas.  An  informal 
report1  lists  these  references  which  is  intended  to  encourage  others  to  become  informed  of  the 
many  past  efforts.  We  make  no  claim  to  being  exhaustive,  or  locating  the  definitive  works,  or 
identifying  the  inventors.  The  preparation  was  mainly  for  a  small  informal  seminar  (held 
November  20,21  1981  at  Biwa-Ko  in  Japan)  where  we  chose  to  list  experiments  and  simulations 
where  the  time-dependent  plasma-wall  behavior  (say,  near  and  played  an  important 
role  in  the  overall  plasma  behavior.  The  work  to  follow  here  has  a  similar  built-in  slant: 


Figure  1.  Model  for  floating  wall  sheath  simulations.  The  RAM  model  has 

particles  only  In  region  I,  0  <x  <L/2,  reflected  at  x  ■  0  and  absorbed  at  x  ■  L/2. 
The  charge  density  Is  mirrored  about  x  ■  L/2  (but  not  Inverted)  with  the  potential 
solved  for  over  x  »  0  to  x  ■  L  (no  particles  In  region  II,  L/2  <x  <L) .  E  ■  0  at 
x  •  0,  L/2,L.  The  density  of  particles  collected  at  x  »  0,  L/2  is  doubled.  The  RAR 
model  has  particles  over  the  whole  region,  x  -  0  to  x  »  L,  with  particles  reflected 
at  x  »  0,  absorbed  at  x  ■  L/2  and  reflected  at  x  ■  L.  Both  RAM  and  RAR  models 
produce  potentials  that  are  symmetric  about  x  *  L/2.  The  model  is  periodic,  period  L. 
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caveat  emptor.  The  references  given  by  Name  (date)  may  be  found  in  the  referencell]. 

(4)  Initial  Objective  -  Sheath  Formation,  and  Steady  State:  A  first  object  is  to  look 
carefully  at  the  formation  of  the  sheath  and  subsequent  steady  state  (if  there  is  one)  in  a 
plasma  near  a  planar  absorbing  wall  or  probe.  The  steady  state  commonly  found  in  time- 
independent  analyses  uses  warm  electrons  which  follow  a  Boltzmann  relation  (implying  mr— 0) 
and  cold  ions.  Solutions  are  found  only  for  ions  drifting  into  the  sheath  regions  with  a  velocity 
equal  to  or  greater  than  the  acoustic  velocity  vs;  ions  are  presumably  accelerated  to  v,  in  a 
pre-sheath  region  many  Debye  lengths  long  called  the  pre-sheath;  there  is  negligible  electron 
density  next  to  the  wall  where  Child’s  Law  applies  (for  space-charge-limited  ion  currents) 
describing  ion  acceleration  into  the  wall.  An  objective  is  to  compare  solutions  of  this  nature 
with  the  time-average  of  our  time-dependent  results. 

(5)  Initial  Model:  The  model  chosen  is  a  uniform  plasma  put  in  contact  with  an  absorb¬ 
ing  wall  at  f  — 0.  Cipolla  and  Silevitch  (1981)  solve  this  time-dependent  model  with  a  fluid-ion 
Boltzmann-electron  analysis.  Denavit  (1979)  treats  a  uniform  warm  plasma  slab  which  is 
allowed  to  expand  freely  after  t  —  0,  using  ion  and  electron  particle  simulation.  These  two  one¬ 
dimensional  models  are  related  in  that  the  absorbing  wall  which  is  unconnected  (floating)  col¬ 
lects  charge  and  creates  an  electric  field  in  the  plasma  similar  to  that  of  the  charge  that  has 
escaped  from  the  expanding  plasma.  Both  articles  provide  thorough  physical  and  mathematical 
descriptions  and  extensive  references  to  prior  work  which,  of  course,  extends  back  50  years  to 
Langmuir. 

The  model  which  we  simulate  is  an  initially  spatially  uniform  warm  plasma  bounded  by  an 
absorbing  wall  (connected  by  a  resistor  R  to  ground,  but  with  R  — floating)  at  one  side  and 
a  charge  reflecting  plane  at  the  other  side,  as  shown  in  Figure  1.  The  model  period  is  L.  twice 
the  separation  between  reflecting  plane  and  absorbing  wall;  either  the  plane  or  the  wall  may  be 
considered  the  center  of  the  whole  period.  The  charge  density  is  mirrored  across  the  absorbing 
wall,  but  not  inverted.  £  — 0  at  x  —  0  (no  surface  charge)  and  at  x-L/2  (by  symmetry,  but 
open  to  question).  Both  ions  and  electrons  are  represented  by  particles,  suiting  with  non¬ 
drifting  Maxwellian  distributions  in  v  at  r- 0.  The  floating  wall  absorbs  particles  crossing  the 
wall  boundary  and  remains  charged  with  them.  After  a  long  time  the  plasma  eventually  decays 
to  the  wall;  there  is  no  reinsertion  of  particles.  With  Te  much  larger  then  T,,  the  wall  very 
quickly  becomes  charged  negatively  (in  about  3  r^),  reflecting  most  electrons,  and  accelerating 
the  ions  that  are  nearby,  subsequently  collecting  electrons  and  ions  in  very  nearly  equal 
numbers.  This  behavior  is  expected  from  elementary  physics.  However,  the  wall  potential  4>w 
is  observed  to  have  large  amplitude  fluctuations  at  the  plasma  frequency  with  the  time  average 
of  6  close  to  but  not  exactly  equal  to  the  floating  potential  expected  from  elementary  time- 
independent  analysis.  The  details  of  the  progress  made  so  far  are  now  given. 

(6)  Typical  Parameters  --  Prototype  for  Warm  Electrons,  Cool  Ions,  Te»T,  >0  :  The 
model  length  L  needed  to  observe  the  physics  is  typically  on  the  order  of  or  larger  than  100 
Debye  lengths,  L  >  lOOX#,.  The  cell  size  Ax  needed  for  numerical  subility  is  something  like  a 
third  of  a  Debye  length  or  smaller,  X^/Ax  >  1/3.  These  choices  indicate  use  of  L  >  300Ax  or 
more.  The  number  of  particles  per  cell,  needed  to  keep  the  particle  (shot)  noise  small  is  about 
5  to  10;  the  number  needed  per  Debye  length  to  keep  the  essence  of  the  plasma  physics  is  also 
about  5  to  10;  hence,  the  toul  number  of  particles  needed  is  on  the  order  of  a  few  thousand. 
These  minimal  parameter  requirements  are  met  easily  and  economically;  observation  of  subtle 
effects  requires  many  more  particles. 

The  oscillation  time  scales  are  the  electron  plasma  period  and  the  much  longer  ion  sound 
wave  and  ion  plasma  periods.  Hence,  the  time  step  is  chosen  to  be  on  the  order  1/15  of  a 
plasma  cycle,  o»Ar  —  0.4.  One  long  time  scale  met  is  the  ion  sound  wave  transit  time  over  the 
system  length,  T-L/v,;  the  other  is  the  ion  thermal  transit  time  T -L/v„.  For  mass  ratio 
»!,/m,-400,  electron  thermal  speed  v„  - 1.0,  vKun<t  s  v„  -Jlnjrn,  —  0.05,  so  that  the  transit 


[1]  C.K.  Birdsall  "PUsma-Sheath-Wall  Time-Dependent  Behavior:  An  Informal  Survey."  ERL  Memo. 
UCB/ERL  M82/51  June  22,  1982. 


r 


Figure  2 
of  100, 


-  11  - 


.  Dispersion  diagram  for  L/R  *  100  (L  ■  lOOA^,  R  ■  10)^) 

T»V 


for  mass  ratio 


time  is  f—  \00\D/{v,eyJmflm,  —  100/w^~  16r^  this  is  about  300r^  or  T— 2000  or  about 
5000  time  steps.  This  is  a  long  time  to  run.  One  way  around  is  to  make  the  mass  ratio  smaller 
(at  mjmt  —  100  the  run  length  drops  by  2).  Another  way  is  to  use  an  implicit  particle  moving 
scheme  with  much  longer  time  steps,  to  be  used  only  after  achieving  assurance  that  short 
wavelength  and  high  frequency  effects  may  be  ignored,  an  assumption  not  yet  gained.  Hence, 
using  a  10/usec/particle/step  computer,  with  104  particles  and  5000  steps,  runs  will  be  500 
seconds  long,  about  8  minutes.  Our  initial  runs  are  for  500  steps,  lasting  less  than  a  minute  on 
the  FACOM  230  Fujitsu  machine  at  IPP;  restart  capability  is  being  worked  out,  to  do  500  steps 


at  a  time.  For  ease  of  reading  results,  we  set 


v„  - 1. 


o>p,  - 1 .0  making  X  &  - 1 .0.  Also  using 


(q/m),-- 1.0  makes  e<b/KTe  -  v,t-)  -<2>;  that  is,  <t>  is  already  in  units  of  e<b/KTf. 

Typically  L  -200  so  that  the  active  region  used  is  L/2  —  100A/>,.  With  512  cells, 
X/v/Ax  =5=  2.5. 

(6)  Dispersion  for  the  Simulations:  The  w  -  k  dispersion  is  shown  in  Figure  2  for 
mjm,  -100,  T,«Tt.  The  lower  branch  is  the  ion  acoustic  wave.  The  upper  branch  is  the 
Gould-Trivelpiece  wave  which  is  forced  on  the  system  by  replacing 


by 


(Or 


)' 


•0»n 


k2 


k*+ky 


k  -  2.405 
R 


Here  R  is  inferred  to  be  the  radius  of  disk  charges  rather  than  sheet  charges  (for  which 
R  — °°).  The  reason  for  forcing  this  dispersion  is  somewhat  obscure,  having  to  do  with 
observing  large  amplitude  oscillations  in  <&»„,/  (to  be  described  shortly)  and  wondering  if  these 
oscillations  drive  large  amplitude  waves,  e.g.,  solitons,  which  require  normal  (nondecay  type) 
dispersion. 

Separate  simulations  with  one  species  with  excitation  localized  in  space  and  lasting  about 
rp  indeed  show  nonlinear  waves,  with  zero  crossings  propagating  away  from  the  region  of  exci¬ 
tation  at  very  nearly 

v-  (up2/kj+  v,e2)''-  =  vcr 

which  is  the  (<o,k— 0)  speed  for  Gould-Trivelpiece  waves. 

Our  simulations  are  done  with  L/R  — 10  (L-  lOOA^,  R  —  lOX^,  )  and  also  with  R  — 
both  "radii"  produce  grossly  the  same  results.  Of  course,  for  /?—<»,  the  upper  branch  is  the 
usual  Langmuir  wave,  with  anomalous  (decay  type)  dispersion. 

(7)  Typical  Results,  for  m,/me  - 100:  Typical  results  are  shown  in  Figure  3  for 
ntf/m,  -100,  and  TJTe  - 10~4.  These  are  snapshots  taken  at  times  50,100,... ,400,  roughly 
every  87^  or  0.87^  (a)  is  <t>  vs  x;  the  slanted  line  drawn  on  the  output  moves  away  from  the 
absorbing  wall  at  speed  v*  —  vteyjmf/m,  —0.1,  at  about  the  first  zero-crossing  of  <J>  (a  little 
meaningless  as  d>  is  not  fixed  at  either  x  -  0orL/2).  (b)  is  the  electron  velocity  distribution 
/(v,)  over  the  whole  plasma  showing  loss  of  the  faster  electrons  and  a  squaring  up  of  /(vf). 
(c)  is  electron  phase  space  showing  both  a  cooling  and  shrinking  near  the  wall,  with  boundary 
region  also  moving  to  the  left  at  about  vf.  (d)  is  /(v,)  with  tick  marks  at  v„  showing  begin¬ 
nings  of  beams  late  in  time,  (e)  is  ion  phase  space,  same  y,  tick  marks,  showing  expansion  of 
the  ion  acceleration  region  at  speed  vs. 

Clearly  it  is  desirable  to  have  these  pictures  developed  over  much  longer  time.  Perhaps 
the  boundaries  moving  at  vs  reach  x-0  (the  center  of  the  whole  period)  and  then  the  ions 
drift  to  the  wall  at  some  speed  on  the  order  of  v5  or  v„.  Runs  with  m,/me  —  1,  vs  - 1.0  to  time 
200  and  with  4,  Vj-0.5,  to  time  400  (both  are  two  ion  sound  transit  times)  show 

something  like  this  but  with  mean  ion  speed  still  increasing  in  space  from  x-0  to  x -L/2. 
However,  as  the  electrons  cool  considerably  (and  hence  vs  decreases)  during  the  run  (no 
reinsertion  or  heating  is  used),  the  ion  sound  may  not  have  yet  reached  x-0.  Very  long  runs 
are  needed. 


(8)  Wall  potential  oscillations:  The  potential  of  the  floating  wall,  is  shown  in  Fig¬ 
ure  4  as  a  function  of  time,  (over  400/ 2 =*647^).  The  dominant  short  period  oscillations 
(period  =20=07^)  are  those  of  mode  2  or  3  at  kk^  <  0.1,  where  no  damping  is  expected. 


POTENTIAL 


V.  vs.  x 


V;  VS.  X 


(d) 


(e) 


Figure  3.  Early  time  history  of  the  formation  of  the  sheath  at  a  floating  wall. 

The  distribution  functions  f(v)  are  averaged  over  all  of  the  particles.  The  slanted 
lines  correspond  to  ion  sound  velocity,  showing  a  sound  wave  propagating  away 
from  the  absorbing  wall  at  roughly  sound  speed,  in  the  potential  and  in  f(v) . 


f 


Figure  4.  Potential  at  the  absorbing  wall  versus  time.  The  fast  oscillations 
are  plasma  oscillations.  Curve  (a)  is  the  kinetic  energy  of  an  ion  absorbed  at 
the  wall  (m.v^/2) .  Curve  (b)  is  a  running  average  of  0wa^ •  “iaioum  value 

near  t  •  60  is  about 


near  t  •  60  is  about 

-^lnCn^/n^)  +  0.25  *  -  2.05  (for  m^/n^  *  100) 
as  expected  from/ theory,  for  T  »  T^.  After  t  ■  100,  electrons  with  v>v£ 
have  free  streamed  at  least  0  to  L,  with  many  collected,  so  that  0  rises. 
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Figure  5.  Model  results  using  Langmuir  dispersion  (&*oo,  Id  sheets)  first  for 
H  ■  N  -  4096,  then  32,768  particles,  mass  ratio  400,  T  Th«  ®«an  0min 

near  t  ■  30  to  80  is  close  to  the  expected  value  from  statics,  -  2.75.  Note 
the  diminution  of  plasma  oscillation  amplitude  as  N  Increases,  but  with  al  most 
no  change  in  p 
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The  points  on  curve  (a)  are  for  the  kinetic  energy  of  the  ions  absorbed  at  the  wall;  the  circles 
are  at  r  — 50,  100,  •  •  •  400,  and  the  x’s  are  averages  over  /  —  0  to  50,  50  to  100,  •  ■  •  350  to 
400.  Curve  (b)  is  an  estimate  for  the  running  time  average  of  <t>.  The  difference  between  (a) 
and  (b)  starts  at  about  b(e&/KTt)  =  1  and  drops  to  about  0.5  later;  this  is  expected  to  be  zero 
in  time-independent  analysis  (the  ions  are  accelerated  from  rest  in  this  time  period).  Due  to 
b<t»/bt  s*  0,  we  cannot  expect  (#mv,*  +  0,^()  to  be  constant. 

In  runs  with  radius  ° °  (Langmuir  dispersion),  the  dominant  oscillations  in  <t>  are  at 

utp  (as  determined  from  measurements  using  different  mass  ratios,  down  to  m,/me  - 1)  and  not 
at  (w^  +  l^v,2)'- ,  hence  at  k\D  ^  0.1  where  no  damping  is  expected.  Typical  results  are  shown 
in  Figure  5,  and  (essentially  the  surface  charge  density),  as  dependent  on  time,  for 
mjm,  -  400,  T,  »  T,.  The  results  for  S,  -  Nt  -4096  show  larger  oscillations  in  than 
seen  for  N,  —  A'’,  —  32,768  particles.  However,  e<<b>/KTfs=2  near  t  —  200  in  both  runs;  over¬ 
laying  the  two  curves  of  shows  only  smoothing  as  .V  is  increased,  with  no  gross  change. 
The  transient  in  shows  electron  collection  until  /=  30,  when  the  first  ion  is  collected,  then 
recovery  to  an  oscillating  state.  We  suspect  that  there  may  be  low  frequencies  (e.g.,  <o  <  up) 
present  m  and  p^u  requiring  longer  runs  to  resolve.  We  suspect  further  that  there  may 
be  intermittent  electron-then-ion  oscillation  bursts  as  postulated  by  us  when  working  with  oscil¬ 
lating  virtual  cathodes  in  electron  diodes  (Birdsall  and  Bridges.  1961).  These  were  seen  by 
Cutler  experimentally  (1964)  and  examined  in  detail  by  Culter  and  Burger  (1966)  in  simula¬ 
tions  and  experiments.  See  also  articles  by  Burger  (1964,  1965,  1967)  and  others  as  given  in 
the  survey  mentioned  earlier.  Resolution  of  these  effects  usnig  simulation  (in  part  repetition  of 
Burger's  pioneering  work,  but  with  improved  techniques  of  simulation,  and  more  detailed  diag¬ 
nostics)  requires  many  more  particles  and  much  longer  runs,  and  introduction  of  a  source  and 
or  re-insertion. 

The  overall  behavior  of  a  drop  in  4  to  a  minimum,  followed  by  a  slow  rise  rather  than 
remaining  at  a  minimum,  is  probably  due  to  the  early  loss  of  the  faster  electrons  (as  noted  by  J. 
M.  Dawson  in  a  seminar).  Hence,  with  a  replenishing  source  supplying  a  full  Maxwellian,  the 
slow  rise  would  not  be  expected.  Thus,  comparison  with  steady  state  expectations  (such  as 
Emmert  et  al  [2])  are  to  be  made  at  rough  checks  with  [2]  are  good. 

(9)  Energy  Histories  —  Ion  heating:  The  active  particle  energy  histories  are  shown  in 

Figure  6  to  accompany  the  earlier  Figures  2,  3,  4.  The  active  particles  are  in  the  region 
0  <  x  <  L/2  (no  particles  L/2  <  x  <  L)  and  the  energies  do  not  include  the  absorbed  particles. 
The  kinetic  and  thermal  energy  of  the  electrons  halves,  as  the  electrons  give  energy  to  the  ions 
and  the  fields  and  are  collected;  at  r»  400  about  16  percent  of  the  electrons  (and  ions)  have 
been  absorbed.  The  drift  energy  of  the  electrons  is  relatively  small,  oscillating  in  time  like 
*„•«.  The  ion  energies  all  increase  in  time  to  about  one  tenth  that  lost  by  the  electrons.  (The 
scales  are  correct  as  ££,(f™0)«  Vsm  n  vw2  —  50.  The  field  energy  is  not  shown,  as  it  was 
obtained  incompletely.  Hence,  there  is  no  total  energy  check  as  yet.)  The  ion  thermal  energy 
shows  an  initial  rise,  apparently  exponential,  from  TJ  T,  -  10-4  to  about  10":.  If  the  growth  in 
T,  is  taken  as  exp(2yt),  then  y  is  measured  to  be  0.77  At  r-0,  kD,/Xx  is  0.0256.  grow¬ 
ing  exponentially  to  0.256  at  TJT,  - 10-2,  and  a  little  larger  later  in  time.  This  initial  heating 
was  observed  further  -hown  in  Figure  7,  for  initial  values  of  T/T,  up  to  about  0.25,  for 
which  -  1"  .  'aps  this  growth  is  nonphysical,  related  to  the  well-known  heating 

due  to  aliasing  obs.r  one-species  Maxwellian  for  kD/±x  <  0.3. 

(10)  Partick  1  •  T  les  are  lost  from  the  plasma  and  absorbed  by  the  wail  at  rates 

shown  in  Figure  8  i  ^  of  Figures  2,  3,  4,  6.  N  is  the  number  of  electrons  or  ions 
remaining  in  the  plasi....  -  measured  at  times  0,  SO,  100,  -  400.  (The  rapid  fluctuations  are 

not  shown.)  The  average  difference  <  /V,  -  A,  >  is  about  15  particles  over  most  of  the  run,  out 
of  2  x  4096  initially,  or  about  0.2  percent.  648  electrons  and  633  ions  are  absorbed,  about  16 
percent  of  each.  The  average  ion  flux  is  thus  A.V/A  T -  633/400  -  1.58.  If  the  ion  flux  is 

12!  O.  A.  Emmsrt,  R.  M.  Wisisnd,  A.  T.  Meus*.  i.  W  Davidson.  "Electnc  shsoth  and  prtalwsih  m  a  coSi- 
donlssa.  Unit*  ion  temperature  pbama*.  Phyt.  /foh.23.  pp.SOJ  - 112.  April  19S0 
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Figure  6.  Various  energies  as  a  function  of  time.  The  loss  of  electrons  (especially 
die  faster  ones)  causes  a  drop  in  KE  by  a  factor  of  2.  The  ion  kinetic  and  thermal 
energies  show  an  early  exponential  rise  (the  latter  is  on  a  log.  scale,  covering 
almost  3  decades)  which  may  be  numerical  in  origin  (self-heating)  due  to  very  small  initial 
Xqj/  x»  «»  a  result  T*/T#  rises  by  about  100  to  where  x  •  0.25. 
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Figure  7.  More  detail*  on  the  initial  rise  of  ion  thermal  energy,  a*  the  Initial 
is  increased,  versus  time.  The  ion  energy  scale  is  logarithmic.  All  for 
■j/®^  •  100,  AjjV  x  ■  2.56.  The  initial  growth  (probably  numerical)  turns  to 
decay  at  largeT^ ,  which  is  probably  mostly  loss  of  ions  (not  coolingX 


taken  as  ( N/L)vt  -  (4096/100X0.1)  -4.096,  the  result  is  too  large,  mostly  due  to  using 
nhVl,  “  W/L  which  is  much  too  large;  the  (low  velocity  at  the  wall  is  1  to  2  times  vs. 

(11)  Comparisons  with  simple  time-independent  analysis:  It  is  useful  to  make  comparis¬ 
ons  of  these  results  with  those  obtained  from  time-independent  analysis.  We  have  already  seen 
one  such,  between  ( ft  m,  v, 2)  and  (q,<l>)  at  the  wall  in  Figure  4.  Let  T,«T,,  essentially  treat¬ 
ing  the  ions  as  cold.  Assume  that  the  ions  are  accelerated  from  rest  up  to  energy  q6vat,  so  that 

C6m,v,2)wll  -  (-q,q>)waii 

not  wholly  inconsistent  with  Figure  3(a)  and  (e).  Let  the  electrons  be  warm,  following  a  sim¬ 
ple  Boltzmann  factor  density  law.  Assume  that  the  electron  and  ion  wall  (luxes  are  the  same, 

r,  -  r„  as 

T,  -  (~)  {vltyf^)txv{e6l KTt)  -  n0v,  -  T, 

This  relation  does  not  account  for  density  at  the  wall  decreasing  as  electrons  are  repelled  or  as 
ions  are  accelerated.  Replacing  <i>  with  (—  ft m,  v,  Vtf, )  and  KT,  with  m,v„2  produces  the  solu¬ 
tion  for  (v,/vlt)wau  vs  mjm,  shown  in  Figure  9  labelled  simple  theory.  The  simulation  results 
are  shown  as  circles.  The  sound  velocity  (v,  —  vtfS[mJm,)  is  also  shown,  substantially  smaller 
than  the  simulations.  The  simulation  results  imply  variation  as  ( mjm,)>  where  p  =  0.4,  not 
square  root.  The  analysis  comes  from  discussions  with  Dr.  Allan  Boozer  at  IPP,  Nagoya. 

Using  the  same  simple  theory,  <f>wan  is  obtained  and  shown  in  Figure  10,  compared  with 
<4>>  from  simulations  as  indicated  in  Figure  S  (late  in  time,  after  some  of  the  faster  electrons 
have  left).  The  time-independent  result  is  about  ft  unit  smaller  than  that  measured;  this 
result  is  similar  to  that  shown  earlier  in  Figure  4  for  mjm,  - 100.  <<t»>  from  simulations  here 
is  measured  qfter  the  initial  transient,  near  the  end  of  the  various  runs  (after  500  or  1000 
steps).  Another  theory  from  [2]  (their  Eq.  (35)  or  their  Figure  4,  for  TJT,  —  0)  is 

■~r  =  ~  ft  in—  +  0.25 
KT  m, 

also  shown  in  Figure  10.  Using  <<&>  as  4^,  before  too  much  loss  of  the  faster  electrons, 
would  lead  to  better  agreement  with  this  latter  theory.  We  have  also  done  (Cal  Tech)  the  dou¬ 
ble  sheath  boundary  (floating  boundaries  at  x  —  0,  L),  not  necessarily  symmetric  about  L/2,  and 
checked  our  early  time  4>mm  vs  mjm,  with  12]  and  found  good  agreement. 

(12)  Ion  wall  speeds  as  TJT,  is  increased:  The  results  with  warm  ions  are  only  partially 
explored.  One  pressing  question  is  the  distribution  of  energies  of  the  collected  ions  for 
T,  =  T„  as  for  fusion  first  walls.  This  is  answered  in  Figure  11.  For  TJT,  ^  0.1,  (after  the 
initial  transient)  all  the  ions  are  collected  at  1  to  2v,,  independent  of  TJT,.  For  larger  TJT„ 
writing  vt-v„Jm,/m,  +  TJT,,  the  mean  value  of  ivjv,,)^  also  is  just  larger  than  v5. 
However,  as  soon  as  TJT,  exceeds  about  1,  the  spread  of  velocities  of  collected  ions  increases. 
The  quality  of  these  runs  decreases  as  TJT,  increases  because  of  the  more  rapid  decrease  of 
particles  remaining  in  the  plasma.  A  source  or  reinsertion  of  particles  is  needed. 
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Flgure  8.  Total  number  of  Iona  and  electrons  remaining  in  the  system  as  a  function 
of  time',  for  4096  of  each  at  t  -  0.  Goes  with  Figures  2, 3,4,6.  Data  taken 
at  points  shown;  later  data  taken  continuously  is  much  the  same  (i.e.,  no  clear- 
cut  oscillations).  (The  data  needed  for  such  is  flux,  like  dN/dt,  not  done  yet.) 


Figure  9.  Speed  of  the  lone  at  absorption  versus  mess  ratio,  with  simulations 
shown  as  open  circles. 


Figure  10.  tfall  potential  taken  at  late  times  (after  the  faster  e's  are  absorbed) 
a  function  of  mass  ratio}  the  theory  from  Ref.  (2)  for  Tg»  ^  is  also  shown. 


Figure  11  Ion  velocity  at  absorption  as  a  function  of  T^/Tc  for  m./m  -  400. 
The  data  become  lass  reliable  for  increasing  T.  due  to  more  rapid  loss  of  ions. 


SECTION  II:  CODE  DEVELOPMENT 


Plasma  Diode:  1-D  Vlasov  Simulation  (GASBAG  Code) 

Wm.  S.  Lawson 


Purpose 

GASBAG  is  designed  to  solve  electrostatic  problems  in  one  dimension  with  non-periodic 
boundary  conditions.  It  solves  the  Vlasov  equation  on  a  grid  vt  -  x.  It  is  being  used  to  provide 
an  independent  check  of  results  obtained  by  particle  methods  which  are  suspected  of  being 
numerical  or  noise  effects. 

Method 

A  non-periodic  Vlasov  simulation  poses  a  major  problem  for  grid  solution,  because  there 
can  be  jump  discontinuities  in  the  distribution  function,  which  no  finite  difference  equation  can 
adequately  handle.  To  solve  this  problem,  a  (hopefully)  continuous  distribution  function  is 
introduced,  and  a  string  of  test  particles  is  used  to  keep  track  of  where  the  discontinuity  lies. 
This  is  a  generalization  of  the  Waterbag  model  used  by  Berk  and  Roberts1.  This  string  of  parti¬ 
cles  can  be  thought  of  as  a  line  dividing  a  real  distribution  function  from  a  fictitious  distribution 
function. 

These  particles  are  moved  using  the  standard  leap-frog  method,  which  is  second  order  in 
time.  A  scheme  for  moving  the  distribution  function  which  is  compatible  with  this  method  of 
moving  the  test  particles  is  the  splitting  scheme  used  by  Cheng  and  Knorr  .  This  splitting 
scheme  has  the  advantage  that  it  can  be  used  with  any  one  dimensional  method  of  interpolation 
(or  finite  difference  approximations  to  the  derivative). 

Current  Status 

The  code  has  been  thoroughly  debugged  for  a  single  species,  and  several  new  diagnostics 
have  been  added.  Work  has  started  on  documentation,  and  a  paper  for  publication.  A  poster 
presentation  of  the  method  was  also  given  at  the  Tenth  Conference  on  the  Simulation  of  Plas¬ 
mas. 

The  code  is  being  used  first  in  a  perturbation  study  of  the  one  dimensional  diode  (i.e.  sin¬ 
gle  species).  The  cathode  emission  distribution  is  perturbed,  and  the  resulting  current  perturba¬ 
tions  are  observed.  This  is  a  first  crude  look  at  how  noise  affects  the  diode.  The  theory  for  this 
is  difficult,  since  the  problem  (even  the  linearized  problem)  is  non-linear. 

Future  Work 

The  GASBAG  code  is  applicable  to  many  interesting  problems.  Some  possibilities  are  the 
study  of  the  filling  in  of  hole  or  trapped  distributions,  and  the  study  of  double  layers. 

References 

1.  H.  L.  Berk  and  K.  V.  Roberts,  Phys.  Fluids  10(1967),  1595. 

K.  V.  Roberts  and  H.  L.  Berk  Phys.  Rev.  Learn  101(1967),  297. 

2.  C.  Z.  Cheng  and  Georg  Knorr,  J.  Comp.  Phys.  22(1976),  330-351. 
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SECTION  II:  CODE  DEVELOPMENT 


B.  Poisson  Solver  with  Boltzmann  Electrons 

V.  A.  Thomas,  A.  B.  Langdon  (Prof.  C.  K.  Birdsall) 


A  nontrivial  problem  in  1-d  simulations  is  to  obtain  4>(x)  from 


d2d>(x) 
dx 2 


-e  [  n,(x)~  n0txp(e<t>(x)/Te)) 


e>0 


(1) 


where  n,(x )  is  known  (from  particle  or  fluid  equations)  and  constants  n0,  Te  are  given.  The 
difficulty  is  with  the  exponential  form  of  the  Boltzmann  electron  density. 


Langdon  suggests  an  iterative  scheme  in  which  we  expand  around  an  approximated  solu¬ 
tion,  as  follows:  Let  5^-*^  —  d>'  where  <b‘  is  known  from  the  last  timestep  or  last  iterate. 
Replace  <t>  with  <f>'  +  8d>  in  Eq.(l)  to  obtain  (expanding  the  exponent  to  0(8<j>)2) 


d2  8<ft 
dx2 


-  ( 


T, 


)exp (e<t>‘/Te)  Btf>  -  - 


~2  +  e(‘n'  ~  n0txp(~)) 


+  0(8d>)2 


(2) 


Note  that  the  left-hand  side  displays  the  local  shielding  with  its  kD~2  multiplier.  The  right-hand 
side  is  treated  as  a  known  source  (ignoring  terms  0(h<t>)2 )  so  that  6<j>  may  be  obtained  directly 
from  Eq.(2).  Then  having  5 <f>,  use  +  repeat  Eq.(2)  and  iterate  this  process 

through  to  convergence.  Successive  iterations  yield  decreasing  8<£’s,  the  corrections  becoming 
ignorable  as  Eq.(l)  is  more  closely  satisfied:  <t>'  —  <t>.  Care  is  needed  particularly  (i)  at  boun¬ 
daries,  (ii)  concerning  overall  neutrality,  and  (iii)  possible  multiple  solutions,  as  Eq.(l)  is  non¬ 
linear.  It  is  worth  noting  that  the  expansion  scheme  used  in  obataining  Eq.(2)  could  be  gen¬ 
eralized  to  any  ne  ( <f> )  i.e.,  is  not  restricted  to  only  exponential  (Boltzmann)  electrons. 

Thomas  uses  a  different  scheme,  taken  from  Mason  [1];  this  second  "direct"  scheme  is 
also  iterative  but  does  not  involve  an  expansion  around  an  approximate  solution.  Start  with 


V£ 


,exp(-^!~)  -  n, 

*  e 


e  >0 


(3) 


In  1-d,  take  dfdtx  of  Eq.(3)  to  obtain 

d2E 
dx2 


“  —el 


e  d6  ,  e  <b ,  a  n 

"•  X  eip(T  ■  ~ 


ed> 


d  n 


(4) 


Replacing  d6/dx  with  and  n0tx p(4r-)  with  ne  -  (e  n,  -  dE!dx)/e ,  we  obtain  an 
equation  in  £  without  exponentials,  which  may  be  written  as 


£ 


d2E 
dx 2 


e  /  dE  \ 
—  Ken,-—) 
Tf  dx 


(5) 


This  equation  is  put  in  finite  difference  form  at  timestep  or  iterate  number  n  on  grids  num¬ 
bered  1  NGy  as 


Ej"  -  G(  £,”,£;". ...£vc".  other  knowns  )  (6) 

where  G  does  nor  contain  £/.  This  means  that  dE/dx  in  Eq.(5)  comes  from  points  j±  1,  and 

It]  R.  Mason,  'Computer  Simulation  of  Ion-Acoustic  Shocks.  The  Diaphragm  Problem". Phvs.  Flwds.14. 
pp.1943-1958.  Sept.  1971. 


-23- 


the  Ej  in  d2E/dx2  is  shifted  to  the  left-hand  side  of  Eq.(6).  Using  previous  iterate  values  on 
the  right-hand  side,  Eq.(6)  with  this  shift  is  then  solved  for  all  Ey,  assume  either  E"  -  0  for 
all  j  or  use  £/~‘  to  start.  Mason  uses  £  -  0  at  the  boundaries  and  "reflection  boundary  condi¬ 
tions"  on  the  density,  i.e.,  at  the  left  edge  j  —  1,  he  sets  nJm0—nJm 2,  and  similarly  for  the  right 
edge. 

Thomas’  problem  has  no  net  electric  field  i.e.,  4 or  equivalently 

f  E-dx  -  0  ; 

hence,  after  predicting  all  £/,  he  uses_(  £/  -£/)  as  his  next  guess:  in  this  way  he  obtains 
convergence  to  Ej  with  zero  average  (£)— 0). 

The  choice  of  which  iterative  scheme  to  use  (the  expansion  scheme  such  as  in  Eq.(2),  or 
the  direct  scheme  of  Eq.(5)  )  is  of  course  dependent  on  the  specific  problem  at  hand.  However 
some  general  statements  can  be  made  about  their  relative  merits.  One  important  consideration, 
for  instance,  might  be  the  scheme’s  convergence  rate  to  a  solution.  For  the  sample  problems 
one  of  us  has  tried  (V.  T.),  the  first  method  of  solution  is  significantly  faster.  This  increased 
speed  does  become  noticeable  even  for  one  dimensional  simulations.  This  conclusion  is  con¬ 
sistent  with  the  conclusions  of  Hockney  and  Eastwood  [2],  who  showed  that  the  first  (expan¬ 
sion)  method  of  solution  has  more  desirable  convergence  properties  than  methods  similar  to 
the  second  (direct)  method  presented  here.  For  higher  dimensional  problems  the  first  method 
also  has  a  natural  extention  and  the  improved  convergence  associated  with  that  method  might 
be  decisive. 


(21  K  W.  Hockney  and  J.  W.  Eastwood,  Computer  Simulation  Using  Particles,  Sec.  6.3  (McGraw-Hitl  Interna¬ 
tional,  I9tl). 
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SECTION  III:  SUMMARY  OF  REPORTS,  TALKS,  VISITORS 


Journal  Articles: 


Douglas  Hamed  "Quasineutral  Hybrid  Simulation  of  Macroscopic  Plasma  Phenomena", 
Jour.  Comp.  Phys,47,  pp.452-462,  1982. 


Reports: 

Y-J  Chen,  W.M.  Nevins,  C.K.  Birdsall,  ‘Stabilization  of  the  Lower  Hybrid  Drift  Instability 
by  Resonant  Electrons",  UCRL-88305  (LLNL)  Oct  28,  1982. 


Talks: 

At  American  Physical  Society,  Div.  of  Plasma  Physics,  Nov  1-5,  1982,  New  Orleans.  LA. 
Abstracts  follow  this  page. 


Visitors: 

Dr.  Siegbert  Kuhn,  from  Professor  F.  Cap’s  group  at  the  Institute  for  Theoretical  Physics, 
University  of  Innsbruck,  arrived  in  December.  He  will  be  a  guest  researcher  in  our  group  for 
his  sabbatical  year(1983).  He  will  be  working  on  formalisms  for  handling  finite,  bounded  sys¬ 
tems  and  aiding  us  in  our  early  efforts  at  simulating  diode-like  systems  and  sheath  phenomena. 
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Tacporal  Development  at  a  Plaama  SheathiOacil- 
Utory  vs.  Steady  State.  Charle*  K.  Birdaall .U.Callf  ■ 
Berkeley  CA  94720*— A  uniform  plasma  of  warm  alactrons 
sad  loss  la  placad  la  contact  with  a  particle  absorb¬ 
ing  wall  which  la  electrically  floating.  Many-particle 
Id  electrostatic  simulation  is  used  to  follow  the  mot¬ 
ion.  There  is  no  source  so  that  the  plasma  eventually 
decays  to  the  wall.  The  floating  wall  quickly  becomes 
charged  negatively;  near  the  wall  the  particles  are 
cooled,  electrons  reflected  and  ions  accelerated,  pro¬ 
ducing  nearly  equal  electron  and  ion  particle  fluxes 
to  the  wall.  During  this  transient  state,  the  potential 
drop  and  the  ion  acceleration  region  (rarefaction  wave) 
propagate  at  nearly  v  .  many  X_  into  the  plasma. The 
plasma-wall  potential  oscillates  ecu)  about  a  mean 
value  <#)■  The  ion  velocity  at  the  wall  falls  between 
vs  sod  /2(q/m)  .<0)‘ .  A  wide  range  of  parameters  is  used 
(mj/m  -  1  to  1600,  T./T  -  10-4  to  64,  L/ A--100, 1000) 
Comparison  is  made  to  a  fluid  model  with  the  same  ini¬ 


tial  conditions  but  using  warm  electrons  with  an  iso¬ 
thermal  Boltzmann  distribution  and  cold  ions**  in 
which  the  ions  were  found  to  enter  the  sheath  at  v  . 
•Work  done  at  I. P.P., Nagoya, Japan  and  at  Calif.  Inft. 
Tech.,  Pasadena  CA;  supported  in  pert  by  US  DOE,  0NR. 
**J.W.Cipolia,M.&Silevltch, J.Plamaa  Phys.25,373.1981. 


ION-ION  TWO-STREAM  IN  THERMAL  BARRIERS: 
VincentThomas  ,  U  .  C . Berkeley .  We  present  stu- 
HTes- 37“ ' the- elecTroTtatTe  ion-ion  two-streasi 
Instability  as  applied  to  thermal-barrier 
cells  of  a  tandem  mirror.  The  basis  of  our 
study  is  a  model  proposed  recently  by 
CohenC 1981  ]  for  the  ion  distribution  function 
in  such  a  cell.  The  results  indicate  that 
this  instability  should  not  present  many  dif¬ 
ficulties  for  parameters  of  interest  in  the 
planned  TMX-U  and  MFTF8  tandem  mirror  experi¬ 
ments  at  LLNL.  These  results  are  quite  dif¬ 
ferent  from  an  earlier  study  which  used  a 
less  accurate  ion  distribution  model  taken 
from  Lontano  ,  PekkerandPozzolK  1980  ] . 

Cohen ,  R.  H. ,  "Axial  Potential  Profiles  in 
Thermal-  Barrier  Cells"  NuclearFu- 
310021.239(  1981  ]. 

Lontano ,  M.,  Pekker,  L.  S.  and  Pozzoli,  R. 
"Inverse-loss  Cone  Instability  in  a  Tandem 
Mirror"  Sov. J  .PlasmaPhys  .6 ■ M32(  1930 1  . 


Ion-Cyclotron  '  Instability .  "NIELS  F.  UTaNI, 
U.C.  Berkeley.  Computer  simulations  are  maae 
of  tne  Aifvdn  ion-cyclotron  (A1C)  instability 
using  two  independent  algorithms.  One  algo¬ 
rithm  advances  tne  vector  potential  on  a  1-a 
grid  parallel  to  B  usi'ng  tne  electron  t  x  tt 
term  in  Ampere's  Law,  wnile  the  other  algo- 
ritnm,  based  on  D.  riarned's  quasineutral 
code,  solves  for  the  electric  field  on  a  2-0 
grid  using  the  electron  momentum  equation. 
Linear  growtn  rates  are  found  to  agree  with 
theory  in  tne  zero  parallel  temperature 
limit.  Saturation  and  suosequent  relaxation 
of  tne  velocity  distribution  occur  via  a 
wave-creaking  aecnanlsm. 

Harned ,  D.  3.  To  appear  in  J.  Comp.  Pnys . 
work  supported  by  DOE. 
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